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1 
NANO-OPTIC FILTER ARRAY BASED 
SENSOR 


The present application is a national phase entry of inter- 
national application PCT/US10/27957, and claims priority to 
U.S. Provisional Patent Application Ser. No. 61/161,892, 
filed Mar. 20, 2009, the disclosure of which is hereby incor- 
porated by reference in its entirety. 


BACKGROUND OF THE INVENTION 


Optical spectroscopy technology has been widely used to 
detect, quantify and analyze the characteristics or concentra- 
tion of a physical, chemical or biological target object such as 
a blood sample. This technology can also be used in other in 
vivo chemometric analyses of chemical components of tis- 
sues or organs in a living organism. À variety of spectral 
techniques involve absorption, transmission, reflection, emis- 
sion, and scattering (including elastic and non-elastic) of 
radiations applied to a target sample. The radiations used span 
over a wide range and include UV, Visual, NIR (Near Infra- 
red), SWIR (Short-Wavelength Infrared), MWIR (Medium- 
Wavelength Infrared), and LWIR (Long-Wavelength Infra- 
red) light. 

Optical spectroscopy is also used for highly accurate color 
measurement of various colored materials. Advanced tech- 
niques are used for clinical quantification of blood glucose, 
dissolved oxygen, dissolved carbon dioxide, urea, lactic acid, 
creatine, bicarbonate, electrolytes, protein, albumin, choles- 
terol, triglycerides, bilirubin, heart rate, breathing rate, hema- 
tocrit, and hemoglobin. 

Optical diagnostics using optical spectroscopy allows for 
the ability to obtain chemical and biological information 
without taking a physical specimen, or the ability to obtain 
information in a non-invasive or non-destructive method 
from a physical specimen. The challenge is that the adoption 
of this technology has been limited due to the size of the 
equipment and the associated cost. Therefore, its application 
has often been limited to centralized labs with scaled testing 
protocols. The opportunity now exists to develop a compact 
and low cost spectrometer. Among those previous efforts to 
miniaturize the spectrometer to expand the application of 
optical spectroscopy into broader uses, the planar waveguide- 
based, grating-based, and Fabry-Perot-based techniques have 
been the major approaches. 

One of the issues encountered when trying to miniaturize 
the spectrometer is the resolution degradation. The resolution 
is usually dominated by the optics, especially by the distance 
from the input slit where the input light comes into the system 
to the detector array such as a photo diode array (PDA). The 
shorter the distances, the higher the resolution degradation 
will be. When filters are used, the number of the filters, and 
the shape or bandwidth (often measured in terms of 
FWHM—Full Width Half Maximum) of each filter dominate 
the degradation. A larger number of filters and a narrower 
FWEM would provide a higher resolution. However, there is 
a certain limitation to how narrow the bandwidths ofthe filters 
can be, especially when these filters are fabricated in an array 
configuration. 

There are some other issues in using these optical spectros- 
copy technologies. For example, Berger et al. (U.S. Pat. No. 
5,615,673) and Yang et al. (U.S. Pat. No. 6,167,290) each 
describe a Raman spectroscopic system designed for trans- 
dermal analysis ofblood components. Xie (U.S. Patent Appli- 
cation Publication No. 2005/0043597) describes a spectral 
analysis system for analyzing blood components using a 
radiation passing through a nail of a finger or toe. In these 
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systems, individual variation in skin or nail properties and in 
blood vessel placement can significantly affect the accuracy 
of the measurements. 


SUMMARY OF THE INVENTION 


Embodiments disclosed herein employ nano-optic devices 
and neural-network-based pattern recognition techniques to 
improve the miniaturized spectral sensing. 

In one embodiment, a device is provided comprising a 
conductive layer including a periodic pattern of elements. The 
elements have shapes and sizes configured such that a trans- 
mittance or reflectance spectrum of the conductive layer has 
a drop at a long-wavelength end. The elements have a period 
configured such that the spectrum has a dip at a Plasmon 
mode resonant wavelength. The spectrum further includes a 
peak between the dip and the drop. 

Methods and the systems to detect, sense or monitor vital 
and health signals using nano-optic filter array based spec- 
trum sensors, and their applications are disclosed. Methods of 
mapping spectrum sensor outputs to health signals or color 
information, and methods and devices used to shape the spec- 
tral responses of the nano-optic filters are also disclosed. 
Structures and configurations of highly conductive materials 
are used to suppress or cutoff the transmission or to enhance 
the reflection. 


BRIEF DESCRIPTION OF THE DRAWINGS 


FIG. 1 is a schematic representation of a spectrum sensor 
using digitized spectral responses of a filter array; 

FIGS. 2A, 2B, and 2C show perspective views of different 
types of plasmonic nano-optic filter devices; 

FIG. 3 is a graphical representation of a processing unit 
within the spectrum sensor; 

FIG. 4 shows the conceptual diagram for the operation of 
the Least Square estimate, equalization, matrices inversion, 
or Moore-Penrose pseudoinverse; 

FIG. 5A is an example of spectral profile of an input light; 

FIG. 5B is an example of a spectral response of a broad 
bandwidth filter; 

FIG. 5C shows an output result without the resolution 
enhancement operation; 

FIG. 5D shows the output with the resolution enhancement 
operation; 

FIG. 6A shows overlaid spectral responses of forty broad- 
band and multi-peak filter sets; 

FIG. 6B shows the individual spectral responses of these 
filters; 

FIG. 7 shows the output from three different LEDs (Red, 
Green and Blue) using the spectrum sensor through resolu- 
tion enhancement operation using digitized spectral response 
data of filters; 

FIG. 8A is a the schematic representation of a wireless 
spectrum sensor; 

FIG. 8B is a prospective representation of a wireless spec- 
trum sensor around a target; 

FIG. 9A illustrates a 2-D array of nano-optic filter structure 
having a rectangular array configuration; 

FIG. 9B illustrates a 2-D array of nano-optic filter structure 
having a hexagonal (or triangular) array configuration; 

FIG. 9C shows an individual slit; 

FIG. 9D illustrates a square array structure having a hori- 
zontal gap that is half the pitch value; 

FIG. 9F-9H show rectangular array structures of various 
vertical and horizontal pitches; 
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FIG. 10 shows low-magnification and high-magnification 
SEM images of various filter structures; 

FIG. 11 shows transmittance spectra of a 1-D line structure 
filter, a 2-D line structure filter, and 2-D array structure filter; 

FIG. 12 shows spectra of 2-D array structures having dif- 
ferent horizontal pitches, illustrating the effects on the trans- 
mittance spectra; 

FIG. 13 shows spectra of 2-D array structures having dif- 
ferent shapes ofapertures illustrating the effect ofthe aperture 
shapes on the spectra; 

FIG. 14 shows the effect of increase in both the horizontal 
and the vertical pitches on the spectra; 

FIG. 15A shows the effect of vertical pitch on the spectra; 

FIG. 15B tabulates the values ofthe pitch and the gaps, and 
the resulting spectral data; 

FIG. 15C shows the effect of the aperture size on the 
spectra; 

FIG. 16 illustrates the effect of the increase in horizontal 
gap on the spectra while the pitch increases proportionally; 

FIG. 17 illustrates the effect of the increase in horizontal 
gap on the spectra while the pitch is kept constant; 

FIG. 18A illustrates the effect of the horizontal pitch on the 
spectra; 

FIG. 18B tabulates the values of the spectra for horizontal 
gaps ranging from 60 nm to 150 nm; 

FIG. 18C illustrates the effect of the horizontal gap on the 
spectra with an increase in horizontal pitch and a fixed verti- 
cal pitch; 

FIG. 19 illustrates transmittance spectra for triangular 
arrays; 

FIG. 20 shows a comparison of transmittance spectra of 
rectangular and triangular array structures; 

FIG. 21 further illustrates advantages of a triangular array 
over a square array; 

FIG. 22 illustrates transmittance spectra from an inverse 
honey comb structure; 

FIG. 23 shows the effect of the aperture area on the spectra; 

FIG. 24 shows measured spectra of various filters; 

FIG. 25 illustrates a system employing the filters disclosed 
herein; and 

FIGS. 26-1 through 26-56 show computer simulations of 
various filter structures. 


DETAILED DESCRIPTION 


Unless otherwise specified, the words “a” or “an” as used 
herein mean “one or more”. The term “light” includes visible 
light as well as UV and IR radiation. 

In FIG. 1, a digital filter spectrum sensor 100 is shown 
containing a set of filters 110, a set or array of detectors 120 
(such as a photo diode array or another suitable photodetector 
array), and a signal processing unit 130. The filters can be 
made of dielectric or metallic materials, can be waveguide 
structures, grating structures, Fabry-Perot etalon structures, 
or plasmonic filter structures. 

The examples of plasmonic filter structures are shown in 
FIGS. 2A, 2B, and 2C. In FIG. 2A, the plasmonic filter device 
210 shows the metallic nanowire or island structure. In FIG. 
2B, the plasmonic filter device 220 shows the metallic film 
with apertures or holes. Examples of such devices 210 and 
220 are described for example in U.S. Patent Application Pub. 
No. 2006/0273245 Al, the disclosure of which is hereby 
incorporated by reference in its entirety. In FIG. 2C, the 
plasmonic filter device 230 shows the metallic embossing 
structures on a metal film. Examples of such devices 230 are 
described for example in U.S. provisional application Ser. 
No. 60/877,660 filed on Dec. 29, 2006, the disclosure of 
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which is hereby incorporated by reference in its entirety. The 
plurality of metal islands, multiple apertures in metal film or 
metal embossing array on a metal film are configured such 
that the incident light is resonant with at least one plasmon 
mode on the structures in said metal film or metal islands. A 
predetermined wavelength will perturb the metallic plas- 
monic structures in surface plasmon energy bands for the 
wavelength selective transmission of light. The filter device is 
located either on a dielectric substrate or in between two 
dielectric layers. 

The spectral response 111 of each filter may show single 
peak as shown in FIG. 5B or multi-peaks as shown in FIGS. 
6A and 6B. The spectral responses of individual filters can be 
measured at the time when the filters are fabricated or during 
the period when the filters are used. The measured spectral 
response data is digitized through sampling in wavelength 
domain over the whole interest wavelength range and quan- 
tization. The sampling rate can be no less than the desired 
wavelength resolution. For example, if the desired resolution 
is 10 nm over the 400 nm range from 380 nm to 780 nm, the 
sampling rate will be no less than 40 times in such a way the 
spectral response of a filter is measured at 485 nm, 495 nm, 
505 nm, ..., and 775 nm, or more. Then the measured data is 
guantized to a digital data at a defined scale. For example, if 
the defined scale is 8-bit, the measured data is converted into 
8-bit digital data. If the defined scale is 16-bit, the measured 
data is converted into 16-bit digital data. The digitized filter 
data now will be stored in the processing unit 130, more 
specifically in the designated memory 331 of the processing 
unit. 

In FIG. 3, a processing unit 300 of the digital filter spec- 
trum sensor 100 is shown as comprising a communication 
unit 340, a processor 310 (such as any suitable logic chip or 
device, such as an ASIC or a computer logic portion), 
memory 320 (such as any suitable volatile or non-volatile 
memory devices or computer components) for program run, 
and designated memory 331 for digitized filter data. The 
communication unit 340 can comprise a serial or parallel 
connector or can be a wireless communication unit such as a 
Bluetooth or Wibree, Zigbee, Wireless USB, or Wi-Fi, or any 
variations of those. The memory 331 is preferably a non- 
volatile memory. In the implementation, this processing unit 
may reside within the same spectrum sensor device 100 or 
may bea separate chip or computer connected to the spectrum 
sensor device 100. The processing operations may be con- 
ducted by software stored in a computer, or by hardware ofthe 
processor 310, depending on the desired configuration or 
whether a computer or a chip is used as a processor. 

In FIG. 1, the input light I,,,,,,,, 140 goes through the filter 
f,.sponse 111 and detector response d,esponse 121 to get an 
output matrix O,,,,,, 122 which is Nx1 array format. The 
detector response d,esponse 121 can be identical to all the 
detectors across the whole detector array. When digitized, the 
filter response f, sponse 111 and detector responsed pouse 121 
are multiplied to combine the spectral responses into a single 
spectral response matrix F such that F,=(d,xf,,, d,xf,, ... 
d af, x), F;7 (di xf5,. d,xf,, ...d,xb,),...,andF,<(d, xf. 
d,xfy, - - dyxfyy). The O,,,,4x now can be expressed as 
FXI,nputi=Ocupus Where n is noise. When the noise n is 
negligible, the unknown input I,,,,,, can be calculated or esti- 
mated through the matrix inversion operation such as 


=p-1 
I, t estimate =F x Output OT 


inpu 


Tput_estimate "F" XO ouput Where F* is psudoinverse of 


Since the pseudoinverse provides a least squares solution 
to a system of linear equations, preferably, pseudoinverse 
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operation may be used for most estimates. FIG. 4 shows how 
the setof digital filter information is used to get an estimate or 
to calculate the input level at each wavelength point. Prefer- 
ably, M (number of estimate point for an input)=N (number of 
filters used) where either matrix inverse or psudoinverse may 
be used. 

In mathematics, the pseudoinverse A* of a matrix A is the 
unique matrix satisfying the following criteria: 

1. AA*A-A; 

2. A"AA"—A" (AT is a weak inverse for the multiplicative 
semigroup); 

3. (AA*)*=A A* (AA? is Hermitian); and 

4. (A*A)*=A*A (ATA is also Hermitian). 

Here M* is the conjugate transpose of a matrix M. For 
matrices whose elements are real numbers instead of complex 
numbers, M*=M7. 

An alternative way to define the pseudoinverse is via a 
limiting process: 


AT =lim(A*A + 61)! A* = limA*(AA* + óD) "1 
6-0 6-0 


These limits exist even if (AA*) and (A*A) do not exist. 
a. Properties 
Pseudoinversion is reversible. It is its own inverse: 
(A*)*=A. 
The pseudoinverse of a zero matrix is its transpose. 
Pseudoinversion commutes with transposition, conjuga- 
tion, and taking the conjugate transpose: 


ADAN, 
Æ*=4", and 


AY AANE 
The pseudoinverse of a scalar multiple of A is the recipro- 
cal multiple of A*: 
(o.4)*-a714* for or. 
If the pseudoinverse of A*A is already known, it may be 
used to compute A": 
A*-(A*4)'A*. 
Likewise, if (AA")" is already known: 
A*-A*(AA*)*. 
b. Special Cases 
Ifthe columns of A are linearly independent, then A*A is 
invertible. In this case, an explicit formula is: 
A*-(A*AY LAH, 


It follows that A" is a left inverse of A: A*A=I. 
If the rows of A are linearly independent, then AA* is 
invertible. In this case, an explicit formula is: 


A -A* (A4*) !.. 
It follows that A" is a right inverse of A: AA*=I. 
If both columns and rows are linearly independent (that is, 


for square nonsingular matrices), the pseudoinverse is just the 
inverse: 


Ara, 


If A and B are such that the product AB is defined and either 
A or B is unitary, then (AB)*=B*A*. If A and B are such that 
the product AB is defined, A is of full column rank, and B is 
of full row rank, then (AB)*=B*A*. The second case here 
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does not cover the first; a unitary matrix must be of full rank, 
but otherwise there is no assumption made on the matrix it 
multiplies. 

It is also possible to define a pseudoinverse for scalars and 
vectors. This amounts to treating these as matrices. The 
pseudoinverse of a scalar x is zero if x is zero and the recip- 
rocal of x otherwise: 


0, 


+3 
P 
xo», 


The pseudoinverse of the null vector is the transposed null 
vector. The pseudoinverse of other vectors is the conjugate 
transposed vector divided by its squared magnitude: 


if x<0; 


otherwise. 


if x=0; 


T 
0, 
xt = x . 
, otherwise. 
xx 


For a proof, simply check that these definitions meet the 
defining criteria for the pseudoinverse. 

c. Finding the Pseudoinverse of a Matrix 

Let k be the rank of a mxn matrix A. Then A can be 
decomposed as A=BC, where B is a mxk-matrix and C is a 
kxn matrix. Then 


A*-C*(CC*)- (B*By !B*. 


If A has full row rank, so that k=m, then B can be chosen to 
be the identity matrix and the formula reduces to A*=A* 
(AA")-!, Similarly, if A has full column rank (that is, kn), 
then A*-(A* A)! A*, 

A computationally simpler way to get the pseudoinverse is 
using the singular value decomposition. 

If A=UZV* is the singular value decomposition of A, then 
A*-VX*U*. For a diagonal matrix such as X, we get the 
pseudoinverse by taking the reciprocal of each non-zero ele- 
ment on the diagonal. 

Optimized approaches exist for calculating the pseudoin- 
verse of block structured matrices. 

Ifa pseudoinverse is already known for a given matrix, and 
the pseudoinverse is desired for a related matrix, the pseudo- 
inverse for the related matrix can be computed using special- 
ized algorithms that may need less work. In particular, if the 
related matrix differs from the original one by only a changed, 
added or deleted row or column, incremental algorithms exist 
that exploit the relationship. 

d. Applications 

The pseudoinverse provides a least squares solution to a 
system of linear equations. 

Given an overdetermined system with independent col- 
umns 


Ax=b, 


we look for a vector x that minimizes 


lláx-bl, 


where || denotes the Euclidean norm. 

The general solution to an inhomogeneous system Ax=b is 
the sum of a particular solution of the inhomogeneous system 
and the general solution of the corresponding homogeneous 
system Ax=0. 
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Lemma: If (AA)? exists, then the solution x can always 
be written as the sum of the pseudoinverse solution of the 
inhomogeneous system and a solution of the homogeneous 
system: 


x= A*(AA*) ! b + (1 — A*(AA* | Ay. 
Proof: 


Ax = AA* (AA* Y | b + Ay — AA" (AA* Ay 
=b+Ay- Ay 
=b. 


Here, the vector y is arbitrary (apart from the dimension- 
ality). In both summands, the pseudoinverse A*(AA)' 
appears. If we write it as A", the eguation looks like this: 


x=A"b+(1-A*A)y. 


The first summand is the pseudoinverse solution. In the 
sense of the least sguares error, it is the best linear approxi- 
mation to the actual solution. This means that the correction 
summand has minimal euclidean norm. The second summand 
represents a solution of the homogeneous system Ax=0, 
because (1-A*A) is the projection on the kernel (null space) 
of A, while (A*A)=A*(AA*Y TA is the projection onto the 
image (range) of A (the space spanned by the column vectors 
of A). The Moore Penrose pseudoinverse is described in more 
detail in Table I which follows the Figures. 

As is usually the case in most signal processing system, 
however, there will be noises in this digital filter spectrum 
sensor system. The noise negatively impacts on estimating 
the input value, reducing the accuracy of the system. To 
remove or reduce the noise effects, simple spatial averaging 
or time averaging can be used. Spatial averaging uses the 
multiple identical sets offilters to receive the input at different 
physical locations on detector array. Time averaging uses 
multiple data reading through the same detector. The multiple 
outputs of the detectors can be averaged, or multiple of the 
input estimates through the matrices inversion can be aver- 
aged. 

Further, when the noise n is not negligible as in most 
practical cases, the unknown input can be estimated through 
various Least Sguare estimate methods with various types of 
constraints, as summarized in the following references: 
Roman Z Morawski, REVIEW ARTICLE, Spectrophotomet- 
ric applications of digital signal processing, Meas. Sci. Tech- 
nol. 17 (2006) R117-R144, and Cheng-Chun Chang and 
Heung-No Lee, On the estimation of target spectrum for filter 
array based spectrometers, 21 Jan. 2008/Vol. 16, No. 2/OP- 
TICS EXPRESS 1056, which are incorporated by reference 
in their entirety. 


In the Chang et al. reference, the Least Sguare estimate ° 


operation is further explained in detail as following. The 
transformation between the target spectrum and the CCD- 
sensor outputs is associated by the matrix equation 

r=Hs+n, (1) 
where the dimensionalities of r, H, s, and n are Nx1, NxM, 
Mx1, and Nx1, respectively. 

Suppose r is an observed signal vector. Working on the 
observation vector, an estimator provides an estimation š of 
the input spectrum by considering all possible source signal- 
vectors s. One criterion which can be used as the starting point 
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is the maximum a posteriori (MAP) rule. The MAP estimator 
is obtained by maximizing the posterior probability, i.e., 


(2) 


SMAP = argmax P(s |r). 


From the Bayes’ rule, the posterior probability can be written 
as P(slr)=P(rls)P(s)/P(r). 

Therefore, there is no information on the source signal such 
that P(s) is uniformly-distributed, the MAP estimator can be 
simplified to the maximum likelihood (ML) estimator. The 
ML estimator maximizes the likelihood function, i.e., 


(3) 


Sui = argmax P(r | s). 
s 


For the filter-array spectrometer, the observed signal vec- 
tor, r, and the source signal vector, s, can be associated by Eq. 
(1) as discussed. Now assume the noise vector n is multivari- 
ate Gaussian with zero mean and covariance matrix R,,, i.e., 
E[n]=0, and E[nn"]<R,,, where the superscript T denotes the 
transpose operation The ML estimator then is obtained by 
maximizing the likelihood function 


exp zs - Hs)" R: (r — Hs)|. © 


P(r|s)= 3 


1 
ONAIR, 


To solve for the estimator, it is equivalent to find the vector s 
which minimizes -2r'R,” 'Hs+s"H"R,, ‘Hs. The solution 
can be found by solving the partial differential eguation 


5 3(-2r7R,,-'Hs+s*H’R,,~'Hs)/ds=0. Providing that the matrix 
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HR, His nonsingular (i.e., inverse exists), the solution of 
this is 
Sa ER HH Ry r. (5) 
Furthermore, if there is no knowledge about the correlation 
of the Gaussian noise vector (or if the elements are mutually 
independent), it is reasonable to substitute the covariant 
matrix R,, by an identity matrix I. Thus the ML estimator, Eq. 
(5), is reduced to the least-squares (LS) estimator, i.e., 
SH Hy HU (6) 
It requires that the inverse of H7H matrix exists. Recall that 
the dimensionality of H is NxM. For solution to exist, M 
needs to be less than or equal to N such that the MxM HH 
matrix is possibly full rank. That is, the number of filters used 
in the filter-array spectrometer needs to be greater than or 
equal to the number of resolved points in the wavelength- 
domain. For the most efficient and practical consideration, 
take M-N, i.e., H 1s a square matrix. Then the LS estimator 


5 can be reduced to 


$, HAY H TEH (7) 

It is worth to mention that, for zero-mean noise, the Ši sr, 
Szs and $,, are unbiased, e.g, E[š,,, (HR, ! H)! H7 
R, 'Hs—s. Therefore, for a fixed unknown source signal vec- 
tor s, one may have the received signal vector r measured 
multiple times over either the temporal or spatial domain. 
This unbiased property ensures the enhancement of estima- 
tion accuracy after averaging operation. The estimation-error 
covariance-matrix of the ML estimator, Eg. (5), can be cal- 
culated and expressed as E[($-s)($8-s) (HFR HYT. We 
note that it is a function of the filter matrix H. Thus, it can 
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show how goodanestimator can be fora particular filter array. 
Although the covariance matrix of system noise R,, is fixed, 
the variance of the estimation error can be amplified by the 
detector sensitivity matrix H. Of interest is the case that H is 
a sguare matrix. Conventionally, the singular value decom- 
position (SVD) is considered as a powerful technigue to deal 
with the noise amplification issue. This method computes the 
inverse of the H matrix based on the singular value decom- 
position where an Eigen value less than a certain threshold 
can be discarded. By exploiting the non-negative nature ofthe 
spectral content, the non-negative constrained least-sguares 
(NNLS) algorithm works particularly well to estimate the 
target spectral. NNLS can be seen as a member of the family 
of the least sguares estimator. NNLS returns the vector S that 
minimizes the norm subject to $70. The original design of the 
algorithm was by C. L. Lawson, and R. J. Hanson. Although 
the NNLS algorithm solves the solution iteratively, the itera- 
tion always converges. 

The input example shown in FIG. 5A is the spectral profile 
of a typical fluorescence lamp. When this input is read or 
measured by the set 40 of the broad bandwidth filters whose 
FWHM is about 60 nm as shown in FIG. 5B, the output 
spectral profile, shown as sguare dot line in FIG. 5C, through 
these 40 filters and detectors is not very close to the original 
input spectral profile. However, when the matrices inversion 
operation is executed, the input estimates, shown as sguare 
dot line in FIG. 5D now become much closer to the original 
input spectral profile. FIG. 6A and FIG. 6B show another 
example of filter set with the broad bandwidth and multi- 
peaks used to detect the spectral profile of LEDs (Red, Green 
and Blue). As shown in FIG. 7, the digital filter spectrum 
sensor with these types of broad bandwidths and multi-peaks 
can rebuild the spectral profile of three LEDs, quite closely to 
the original LEDs’ spectral profile. Noticeably, the band- 
width of the filter used is close to 100 nm, but the digital filter 
spectrum sensor system reconstructs the input signals at 
higher resolution such as 10 nm. Also the known spectral 
information of the LEDs, input lights, can be utilized to 
estimate further closer to the original signals. This color mea- 
surement capability will find lots of useful application areas 
in color measurement or monitoring for flat panel display, 
printed materials, or even in bio applications. The measured 
spectral information can be mapped into color domain data in 
different format such as CIE systems, such as CIELUV or 
CIELAB or CIE 1931 or CIE 1976. 

Because of its compact size and reasonable resolution, the 
digital spectrum sensor can be used in many applications 
where the small size is desired. Such applications may 
include mobile or wearable devices to monitor personal 
health, high resolution color monitoring for the color input, 
display and output devices, and environmental monitoring 
sensors such as water or air quality sensors. For some of those 
applications, wireless solution may be more desirable. As 
shown in FIGS. 8A and 8B, the wireless spectrum sensor has 
one or more embedded wireless units 702 mounted on a 
holder (such as a substrate or a band) along with spectrum 
sensor 704, one or more light sources 701, and a power unit 
703, such as a battery unit or power generator. The holder may 
have embedded circuits connecting these embedded units. 
This wireless digital spectrum sensor may be located adjacent 
to or wrap around the target object 710, such as a human or 
animal body part, an inanimate object, the container of the 
target object to be measured, or the passage or pipeline where 
the target objects (such as any suitable fluid) to be measured 
pass through. If the power unit 703 is a power generator, then 
this power generator may utilize wireless RF power to gen- 
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erate the electrical power and/or it may utilize the flow of said 
object, such as water or other fluid, to generate power. 

FIG. 9A illustrates a 2-D array of nano-optic filter structure 
910 comprising an aluminum layer on a quartz substrate. The 
Al layer has a rectangular array of apertures. The filter struc- 
ture 910 can be fabricated, for example, by electron beam 
photolithography and Al etching process on the quartz sub- 
strate. The filter structure 910 can have dimensions of, for 
example, 50x50 um. FIG. 9B illustrates a 2-D array of nano- 
optic filter structure 920 having a hexagonal (or triangular) 
aperture array configuration. Other array configurations can 
also be employed. 

In the examples shown in FIGS. 9A and 9B, the array 
elements are the slits or apertures in a conductive layer. In 
some other embodiments, the array elements can include 
reflective particles or protrusions. The shapes, sizes, gaps 
between the array elements, and the patterns can be engi- 
neered to affect the transmittance (if the filter is configured as 
a transmissive filter) or reflectance (if the filter is configured 
as a reflector or a reflective filter). 

FIG. 9C shows an individual slit 930, which can be char- 
acterized by a vertical slit width (V-Gap, or y) anda horizontal 
slit length (H-Gap, or x). 

FIG. 9D illustrates a square array structure 940a having an 
H-Gap=0.5xPitch. FIG. 9E illustrates a triangular array struc- 
ture 9405 having an H-Gap=0.5xPitch. The height of the 
triangle is V3/2xPitch. 

FIG. 9F-9H show rectangular array structures 9504, 9505, 
950c of various vertical and horizontal pitches. 

FIG. 10 shows low-magnification and high-magnification 
SEM images of various filter structures of the embodiments 
disclosed herein. 

FIG. 11 shows transmittance spectra 1110, 1120, and 1130, 
respectively of a 1-D line structure filter 1111, a 2-D line 
structure filter 1121, and 2-D array structure filter 1131. From 
these spectral curves 1110, 1120, and 1130 it can be seen that 
the Ay. rear (~530 nm) and X, py (~630 nm) exist at about the 
same wavelength locations regardless of whether the filter 
structure is 1-D or 2-D. The peak intensities of the À paar of 
the 1-D line structure (~0.9) and the 2-D line structure (~1.75) 
are different. The longer wavelength noise of 2-D line struc- 
ture (~5.6) is significantly higher than that of the 1-D line 
structure (~2.2). 

The 2-D array filter structure 1131 as shown has a trans- 
mittance spectrum with a drop at the long wavelength end. 
The drop or cut off can be characterized by a “cut-off” wave- 
length, for example, defined by the wavelength at which the 
transmittance drops to about 50% of its peak value. The drop 
in the spectrum 1130 in conjunction with the dip at A, pyp 
forms another peak at X, py (-750 nm in this case). This 
advantageously narrows the shape of the spectrum 1130 at the 
wavelength region of interest thereby improving the resolu- 
tion of the transmittance spectrum. 

The shape of the transmittance spectrum can be engi- 
neered, for example, by varying the length (x) and width (y) 
of the slits, and the horizontal pitch (X) and vertical pitch (Y). 
X, Y, x, and y are generally comparable to the wavelength 
range of interest (e.g., IR radiation, UV radiation or visible 
light (400-700 nm)). Preferably x and y are smaller than the 
wavelength of interest, and thus the slits are often referred to 
as having “sub-wavelength” dimensions. In a preferred 
embodiment, y<x, and in the filter structure 1131 as shown 
x/X-0.5. For example, X=100 to 1000 nm, x=50 to 500 nm, 
Y=200 to 1000 nm and y=20 to 200 nm. 

In one embodiment, the slits are formed in a highly con- 
ductive layer. The highly conductive layer can comprise, for 
example, a highly conductive metal or an alloy of metals (i.e., 


US 9,395,473 B2 


11 


metal alloy), a highly doped semiconductor, a layer compris- 
ing carbon nanotubes or graphene, or coated with highly 
conductive materials. In some embodiments, the highly con- 
ductive layer can be a combination of multiple layers. 

The highly conductive layer can be disposed between at 
least two dielectric layers (e.g., guartz, silicon dioxide, silicon 
nitride, alumina, etc.), as shown for example in FIG. 26-1. 
The highly conductive layer can have one or more periodic 
pattern of sub-wavelength apertures. When two or more peri- 
odic patterns are included, these periodic patterns can inter- 
lace with each other. Alternatively, the array elements can be 
made of sub-wavelength sized conductive particles or protru- 
sions. 

The conductive layer can be made from a same material 
across the device, although different structures can be dis- 
posed over the conductive layer. Alternatively, at least two 
different highly conductive materials can be used across the 
device. Similarly, the same, or at least two different dielectric 
materials, can be used for the top or bottom dielectric layer. 

The shapes and sizes ofthe apertures or gaps are configured 
such that the transmittance or reflectance spectrum has a drop 
at the long wavelength end. The drop is often referred to as a 
“cut-off” that is characterized by a cut-off wavelength. The 
periods of the array elements can be configured such that the 
transmittance spectrum is minimized or suppressed (having 
dips in the spectral curves) at specified wavelengths. These 
wavelengths can be those resonant with at least one Plasmon 
mode of the periodic patterns. 

As a result of engineering the array structures described 
above, at least one peak in the transmittance or reflectance can 
be located between the dip and the cut off wavelength, as 
shown in the spectrum 1130 in FIG. 11. 

In some embodiments, the thickness of the highly conduc- 
tive layer is thicker than a skin depth, i.e., the light penetration 
depth, of the highly conductive material. 

In some embodiments, the apertures or gaps can have 
sharpened edges to enhance the transmission or reflection of 
light at the peak wavelength. The enhancement of the light 
transmission is realized through, for example, highly con- 
densed charges or fields of the sharpened edges. 

A coating comprising a magnetic material can be included 
to modify the electromagnetic field adjacent the conductive 
layer to thereby affect a Plasmon mode. 

As described in detail below, the narrow widths of the 
apertures or gaps can also enhance the sharpness of dips and 
peaks, thereby reducing the FWHM and improving the reso- 
lution of the filter or reflector. 

In some embodiments, the conductive layer can have at 
least two alternating periods of the periodic pattern, and/or at 
least two alternating shapes or sizes of the array elements. 

In a preferred embodiment, the dielectric layers on both 
sides of the conductive layer are composed of the same mate- 
rial. In addition, dielectric material filling the apertures or 
gaps can also be the same material as the dielectric layers 
sandwiching the conductive layer. By reducing the types of 
interfaces, the number of the Plasmon modes can be reduced. 
Thus, the number of dips and peaks in the transmittance or 
reflectance spectra can be reduced, and the transmittance or 
reflectance can be increased. 

The shapes and sizes of the apertures or gaps on the top 
surface and the bottom surface of the conductive layer can be 
the same, or can be different. 

In some embodiments, a plurality of optical filters dis- 
closed above can be stacked together, wherein a separation 
between the two conductive layers is in the sub-wavelength 
range, or can be larger than a predetermined wavelength. 
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In one embodiment, a method is provided to shape a spec- 
trum. A highly conductive layer can be provided between at 
least two dielectric layers. The highly conductive layer can 
have at least one periodic pattern of sub-wavelength aper- 
tures, or sub-wavelength sized conductive particles. The 
shapes and sizes of the apertures, gaps, or the particles are 
configured such that the transmission of incident light is 
suppressed at the long wavelength portion of the incident 
light. Periods ofthe apertures or particles are configured such 
that the transmission of the incident light is minimized or 
suppressed (having dips in the spectrum) for the wavelengths 
resonant with at least one Plasmon mode of the periodic 
patterns. As a result, at least one peak wavelength of the 
transmission is located between the dip and a cut off wave- 
length that characterizes the drop at the long wavelength 
portion. 

An array of optical filters can be provided, for example, 
including at least two optical filters, or at least four optical 
filters. 

A digital spectrum sensor can be provided. The sensor can 
include a photodetector array, and an array of optical filters. A 
processing unit can be included for estimating input spectral 
resolution enhancement. Each filter of the filter array can be 
optically coupled to a photodetector or a group of photode- 
tectors of the photodetector array. A vector value of the out- 
puts from the photodetectors associated with the respective 
filters is used to analyze properties of a target object or to 
monitor changes in the target object. Methods such as linear 
ormulti-linear or nonlinear estimation, or trained mapping, or 
pattern recognition can be employed in the analysis. 

The spectral response data of the optical filters or com- 
bined spectral responses of each filter and its associated pho- 
todetector can be sampled, guantized and digitized and used 
to analyze properties of the target object, or to estimate of the 
spectral profile of the incident light. 

The vector value can be an array of the outputs from the 
photodetectors and their respectively associated filters. The 
vector value can be a matrix of cross relationships among the 
outputs, a matrix of ratios among the photodetectors and 
filters. The vector value can be a matrix of first or second 
derivatives of the outputs. 

The linear or nonlinear estimation methods can use at least 
one ofthe following operations: matrices inversion, egualiza- 
tion, Moore-Penrose pseudoinverse, least sguare, linear or 
nonlinear regression, neural network, or multilayer neural 
network. 

The photodetector array can use at least two different sizes 
of pixels or at least two different shutter times to enhance a 
dynamic range or to normalize the outputs of the photodetec- 
tor array. Pixel binning method may be used to normalize the 
outputs of the photodetector array. 

The photodetector array can be only partially covered by 
the optical filter array to obtain image data in addition to the 
spectral data, e.g., some of the photodetectors can receive 
unfiltered light. Thus, the device can be configured as a spec- 
tro-photometric sensor conducting spectroscopy and imaging 
at the same time. 

Spatial, temporal, or moving-time averaging processes can 
be used to improve the Signal-to-Noise ratio. The spatial 
averaging includes repeating identical filters over the mul- 
tiple photodetectors and averaging the outputs. The time aver- 
aging process includes reading the output of a photodetector 
repeatedly over the multiple frames and averaging the out- 
puts. 

The digitized spectral response data can be adaptively or 
selectively used for different input lights. The data can be 
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measured and stored in memory as fixed data. The data can 
also be adaptively modified or updated through calibration. 

Input constraints can be used to process the vector value to 
analyze properties of the target object, monitor changes of the 
target object, or estimate a spectral profile of input light. 
Linear or nonlinear estimation, trained mapping, and pattern 
recognition are some of the methods that can be employed in 
the analysis. The known input constraint comprises positivity 
or negativity of the input signals, known spectral information 
of light sources used, boundary conditions of input signals, or 
temporal, spatial, or frequency modulated information. 

A database of pre-tested or pre-measured relationships 
between the vector values and the properties of the target 
object can be stored and used to map or convert the spectral 
response data into the space domain property of the target 
object. 

A digital spectrum sensor can be provided including a 
photodetector array and an array of optical filters. A process- 
ing unit can be used for estimation of inputs or input spectral 
resolution enhancement. The optical filters can include 
Fabry-Perot filters, linear variable filters, thin film filters, 
photonic crystal filters, or nanostructured thin film filters. 
Each filter of the filter array can be optically coupled to a 
photodetector or a group of photodetectors of said photode- 
tector array. 

A digital color sensor can be provided and can include the 
digital spectrum sensor. The vector value of the outputs can be 
converted or mapped into human's three-color sense values or 
tristimulus values. The human’s three-color values are repre- 
sented in the form of International Commission on Illumina- 
tion (CIE) XYZ values, or International Commission on Illu- 
mination (CIE) xy chromaticity values, or xyY, or CIE RGB 
values. 

The stored digitized spectral response data can be used to 
provide the estimated spectral profile of the input light and the 
color values. 

The systems and methods can be used for, e.g., color sens- 
ing for TV, color sensing for ambient light, color sensing for 
printer, or color sensing for skin, food, process monitoring. 

An indirect biochemical material detection apparatus can 
be provided comprising a color sensor, at least one light 
source, a processing unit with memory, wireless or wired 
communication means, at least one color changing material, 
and a light path. The color changing agent changes color 
when exposed to a predetermined target biochemical mate- 
rial. The degree of the color change can be proportional to the 
amount of the biochemical material. 

A method can be provided to indirectly detect biochemical 
material using the apparatus. The method includes exposing 
the color changing agent to a target biochemical object, or to 
the environment where the target biochemical object exists, 
shining light onto the color changing agent, sensing the color 
changes by reading the color sensor outputs, determining the 
amount of contacted target biochemical materials, and trans- 
mitting the information. 

Sensing the color changes can be achieved by filtering a 
reflected, scattered, or transmitted light from the color chang- 
ing agent using the filters or reflectors disclosed herein. The 
filter suppresses a spectrum of the light from the color chang- 
ing agent at a short wavelength region using a Plasmon reso- 
nant mode of a periodic pattern in a conductive layer of the 
color sensor, and suppresses the spectrum at a long wave- 
length region by selecting shapes and sizes of elements in the 
pattern to thereby form a peak in the spectrum. 

A noninvasive health monitoring device can be provided 
comprising at least one digital spectrum sensor, at least one 
light source, a processing unit, a wireless or wired commu- 
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nication unit, optical means for the light path, mechanical 
means for attaching device to part of a body. The light sources 
emit light onto the body, and the sensor detects the reflected, 
scattered, or transmitted light from the body. Vector values of 
outputs of the sensors are converted, mapped, or interpreted 
as the person's vital or health information. The vector values 
can also be transmitted to an information collecting unit. The 
device can be wearable on the body using, for example, a 
band, a clip, or a hook. 

The part of body can be preferably an earlobe, a finger tip, 
a wrist, or an upper arm. A fluorescence material can be 
embedded in said part of body for enhancing the target object 
detection. 

The light from the light source can be modulated by fre- 
guency modulation, amplitude modulation, or code modula- 
tion. The light source outputs can be wavelength multiplexed 
from at least two light emitters. 

The mapping or conversion of vector value to vital or 
health information can be self-trained or self-calibrated for an 
initial use, or after power off-on of said device. A reference 
light source and a reference light paths can used for self 
calibration of the device. 

A multispectral imager or a hyperspectral imager can be 
provided including a photodetector array and an array of 
optical filters. Each of the optical filters can be associated 
with a pixel of the photodetector array. A mosaic pattern of a 
group of the optical filters can be repeated spatially over the 
entire area of the photodetector array. The number of different 
filters within each group can be, for example, at least four, or 
at least twelve. The different spectral outputs of pixels within 
the group can used to define, process, or approximate the 
spectrum of the group and neighboring groups. 

A tunable optical filter or reflector can be provided utilizing 
the structures disclosed herein. A voltage or current source 
can be used to apply a voltage or current to the conductive 
layer and/or dielectric layers. An intensity or wavelength of 
light transmitted through or reflected from the tunable device 
is modulated by the voltage or current. 

The tunable device can be used in, for example, a trans- 
missive, a reflective, a transreflective flat panel display, which 
includes an array of the tunable optical filters. The tunable 
device can also be used in a hyperspectral imager, a photode- 
tector array, an optical modulator, or a Fabry-Perot filter 
comprising two stacked filters in which the distance between 
the two filters is smaller than a predetermined wavelength of 
incident light. The two stacked filters can have different shape 
and configurations. 

In the tunable Fabry-Perot filter, the dielectric materials 
between the two filters can have sub-wavelength structures to 
change the effective refractive index. 

FIG. 12 shows spectra of 2-D array structures 1211, 1221, 
1231 having different horizontal pitches, illustrating the 
effects on the transmittance spectra 1210, 1220, 1230. By 
increasing the horizontal pitch, from example, from 1P to 2P, 
while the horizontal gap is kept at half of the horizontal pitch, 
and the vertical pitch is fixed at about 390 nm and the gap is 
fixed at about 40 nm, it can be shown that the transmittance 
spectra 1210, 1220, 1230 have increasingly higher peaks. The 
slopes of the peaks are different as well. 

FIG. 13 shows spectra of 2-D array structures having dif- 
ferent shapes of apertures illustrating the effect of the aperture 
shapes on the spectra. A simple, rectangular pattern 1311 
comprising intersecting lines, similar to that of the structure 
1121 shown in FIG. 11 and, results in a spectrum 1311 with- 
out a drop at the long wavelength end. A rectangular pattern 
1321, in which the apertures are shaped like crosses with their 
center portions blocked, results in a spectrum 1320 with a 
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drop at the long wavelength end and a relatively broad peak. 
A rectangular pattern 1331, in which the apertures are rhom- 
bus shaped, results in a spectrum 1330 having a desirable 
narrow peak between the dip and the cutoff wavelength (i.e., 
drop). 

FIG. 14 shows the effect of horizontal pitch on the spectra. 
As the horizontal pitch increases from about 200P to about 
300P and then 400P while the vertical pitch and the horizontal 
gap of the slits increase proportionally, the peak locations 
shift red-ward (to longer wavelength) from about 420 nm, to 
about 590 nm, and then to about 760 nm. The peak wave- 
length has a linear relation with horizontal pitch or the diago- 
nal pitch. 

FIG. 15A shows the effect of vertical pitch on the spectra. 
As the vertical pitch increases, the dip shifts towards the long 
wavelength end, and the slopes at the left sides of the peaks 
also change. The pitch is kept the same in horizontal and 
vertical directions, and the horizontal gap remains at 50% of 
the pitch. FIG. 15B tabulates the values of the pitch and the 
gaps, and the resulting spectral data. FIG. 15C shows the 
effect of the aperture size on the spectra. As the horizontal gap 
increases together with the pitch values while the vertical gap 
is kept const, the cut-off wavelengths increase, and the slopes 
at the right sides of the peaks also change. The pitch is kept the 
same in horizontal and vertical directions, and the horizontal 
gap remains at 50% of the pitch. 

FIG. 16 illustrates the effect of the increase in horizontal 
gap on the spectra while the pitch increases proportionally. As 
the vertical gap increases, the transmittance at the right side of 
the peak increases. The noise in the long wavelength range 
can also increase. The dip at about 800 nm may be related to 
the metal (Al) material property. For example, the permittiv- 
ity curve of Al indicates that there may be a high loss in the 
500~800 nm range. 

FIG. 17 illustrates the effect of the increase in horizontal 
gap on the spectra while the pitch is kept constant. The peak 
wavelength increases while the peak broadens slightly. The 
transmitted light per unit aperture (gap) area=Transmittance/ 
Total area of gap, and the total area of gap=(Gap area/unit 
cell)x(# of cell in each filter). 

FIG. 184A illustrates the effect of the horizontal pitch on the 
spectra. While the horizontal pitch increases from 1P to 2P, 
the horizontal gap increases proportionally and is kept at 
0.25xthe horizontal pitch. The vertical pitch is fixed at 390 
nm, and the vertical gap is fixed at 40 nm. The peak of the 
spectrum shifts toward the long wavelength end and becomes 
broader. An optimal configuration may be present, and can be 
determined by simulating a large number of different con- 
figurations. For example, in FIG. 18A it is shown that when 
the horizontal pitch is 585 nm, the peak is substantially nar- 
rower than that for a horizontal pitch of 780 nm while having 
a desirable intensity. 

FIG. 18B tabulates the values of the spectra for horizontal 
gaps ranging from 60 nm to 150 nm. 

FIG. 18C illustrates the effect of the horizontal gap on the 
spectra with an increase in horizontal pitch and a fixed verti- 
cal pitch. The area of the gap is kept at a constant value, with 
a horizontal gap of 98 nm. The transmittance is normalized by 
the hole area. 

FIG. 19 illustrates transmittance spectra for triangular 
arrays. If the area of the aperture is large, the filter acts like a 
polarizer filter, and the intensity slope increases. If the aper- 
ture is small, then the wavelength is limited by the hole size, 
and cut-off mechanism takes place, and thus the intensity 
slope decreases. 

FIG. 20 shows a comparison of transmittance spectra of 
rectangular and triangular array structures. In terms of slit 
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array structure, the hexagonal (or triangular) array shifts the 
peak wavelength for about 60-80 nm toward shorter wave- 
length compared with the rectangular array at the same foot- 
print. 

FIG. 21 further illustrates advantages of a triangular array 
shape over a square array shape. 0,,<0,,; indicates that trian- 
gular arrays can have a sharper peak in its spectrum than 
sguare arrays, resulting in a decreased FWHM. The dip loca- 
tion of the triangular array is shifted for about 13% toward 
shorter wavelength, resulting from that the height of the tri- 
angle is about V3/2=0.87 that of the square. Thus, A, can be 
shifted toward shorter wavelengths at the same footprint by 
engineering the shape of the periodic pattern. 

In the transmittance spectra, the first dip results from a 
Plasmon mode at an aluminum-guartz interface, and the sec- 
ond dip results from a Plasmon mode at an aluminum-air 
interface. Thus, by engineering the refractive indices, for 
example, by selecting the dielectric materials on one side or 
both sides of the conductive layer, and/or selecting the con- 
ductive material itself, the spectral shape can be modified. In 
particular, increasing the types of interfaces can increase the 
number of peaks and dips. Iffewer peaks and dips are desired, 
a single type of dielectric material is desired. This provides 
another design freedom in addition to designing the periodic 
patterns and the shapes and sizes of the array elements. 

FIG. 22 illustrates transmittance spectra from an inverse 
honey comb structure. 

FIG. 23 shows the effect of the aperture area on the spectra. 
As shown, when the slit open area is 59% of the total area of 
the conductive layer, the spectrum does not show a drop in the 
long wavelength end, and the filter can function as a long-pass 
filter or polarizer. When the slit open area is 43% of the total 
area of the conductive layer, the drop or cut off is present. 

FIG. 24 shows measured spectra of various filters. 

FIG. 25 illustrates a system employing the filters disclosed 
herein. The system is configured as a wearable personal 
health monitor, including a strap to attach the system to a 
wrist. Alternatively, with a different mechanical attachment 
device, the system can be attached to a fingertip or an earlobe 
of a body. The system can include a light source that emits 
light onto the skin. They reflected, scattered, or transmitted 
light can be analyzed by a spectrometer on a chip, where the 
spectrometer includes the optical filters disclosed herein. The 
measured raw data can be processed by a data/signal process- 
ing unit. The processing can include a preprocessing, and 
feature extraction and classification. These data can be inter- 
preted as health data, which is transmitted to a wireless mod- 
ule, which in turn sends the health data through wireless (e.g., 
Bluetooth) or wired (e.g., USB) communication device to a 
PC, a personal health monitoring system, or to a cell phone. 

FIGS. 26-1 through 26-56 show computer simulations of 
various device structures. Such simulations can be used to 
design optimized filter configurations. In particular, FIGS. 
26-1 through 26-11 show simulated spectra of various multi- 
layer structures. As illustrated in the cross-sectional views of 
the devices, the distances between the two layers, the shift 
between the two periodic patterns, the different pitches in the 
two layers, etc, can all affect the shape of the spectra signifi- 
cantly. FIGS. 26-12 through 26-24 show the effect of multiple 
periods in the pattern. FIGS. 26-25 through 26-32 show the 
spectra when the periodic pattern is configured as a 1D polar- 
izer or long-pass filter. FIGS. 26-33 through 26-38 show the 
effect of variations in the aperture shapes. 

Although the foregoing refers to particular preferred 
embodiments, it will be understood that the present invention 
is not so limited. It will occur to those of ordinary skill in the 
art that various modifications may be made to the disclosed 
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embodiments and that such modifications are intended to be 

within the scope of the present invention. 
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What is claimed is: 
1. A method of characterizing an incident radiation 
employing filtering, the method comprising: 

providing a spectrum sensor that includes a plurality of 
combinations of a filter device and a photodetector opti- 
cally coupled to the filter device, wherein each filter 
device has a transmittance spectrum that is different 
from transmittance spectra of other filter devices, 
wherein at least one of the filter devices comprises a 
contiguous conductive layer in which the corresponding 
array of periodic patterns are embodied as apertures 
entirely laterally surrounded by the contiguous conduc- 
tive layer, and wherein at least one of the filter devices 
comprises two dielectric layers, the contiguous conduc- 
tive layer is disposed between the two dielectric layers, 
the continuous conductive layer includes a plurality of 
apertures and gaps, and the apertures and gaps are filled 
with a dielectric material that is the same as the two 
dielectric layers to thereby reduce numbers of dips and 
peaks in the spectrum and increase transmittance or 
reflectance; 

providing an incident radiation to each filter device of the 
spectrum sensor; 

detecting, at each photodetector, a transmitted radiation 
from a respective filter device; 

generating measured outputs from the photodetector, 
wherein different transmittance spectra of the filter 
devices provide different weighting to a spectrum of the 
incident radiation; and 
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generating an estimation ofa spectral profile ofthe incident 
radiation from the measured outputs of the photodetec- 
tor employing an estimation method that utilizes the 
different weighting corresponding to the filter devices. 

2. The method of claim 1, further comprising analyzing a 
biochemical material by: 

illuminating a portion of a biochemical material with an 

incident radiation, wherein the biological material is ina 
state that emits a spectrum that is different from the 
incident radiation; 
sensing a color or spectral change in an emitted spectrum 
form the biochemical material with respect to the inci- 
dent radiation using the spectrum sensor; and 

determining an amount of the biochemical material based 
on the color or spectral change. 

3. The method of claim 2, further comprising exposing a 
color or spectral changing agent to the biochemical material 
prior to illumination of the biological material with the inci- 
dent radiation, wherein the color or spectral changing agent 
provides the color or spectral change in the emitted spectrum. 

4. The method of claim 1, wherein the filter devices are 
configured to be simultaneously exposed to the incident 
radiation. 

5. The method of claim 1, wherein each filter device 
includes an array of periodic patterns having a plasmon mode 
resonant wavelength that differs from plasmon mode reso- 
nant wavelengths of other filter devices. 

6. The method ofclaim 1, wherein the transmitted radiation 
has a transmittance spectrum that has a cutoff at a long- 
wavelength end of a spectrum, a dip at a plasmon mode 
resonant wavelength, and a peak between the dip and the 
cutoff. 

7. The method of claim 6, wherein the different dips or 
peaks in the transmittance spectra provide different wave- 
lengths of minimum or maximum weighting for the incident 
radiation to the measured outputs. 

8. The method of claim 7, wherein the estimation method 
utilizes the different wavelengths of minimum or maximum 
weighting corresponding to the different dips or peaks. 
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9. The method of claim 1, further comprising: 

generating an observed signal vector from the measured 

outputs; and 

multiplying a spectral response matrix and the observed 

signal vector to generate a vector representing the esti- 
mation of the spectral profile. 

10. The method of claim 1, wherein the contiguous con- 
ductive layer comprises at least one of a conductive metal, an 
alloy of metals, a highly doped semiconductor, a carbon 
nanotube, a graphene, a material coated with highly conduc- 
tive materials, and multiple layers thereof. 

11. The method of claim 1, wherein the array of periodic 
patterns is a two dimensional array that is repeated along two 
different directions. 

12. The method of claim 1, wherein shapes of the plurality 
of apertures include at least one of a rectangle, a ninety- 
degree zigzag shaped slit, and a cross. 

13. The method of claim 1, wherein the spectrum sensor is 
configured to perform at least one of the following: 

monitoring an ambient light; 

monitoring proximity; and 

color or spectral detection. 

14. The method of claim 1, wherein the spectrum sensor 
further comprises a processing unit configured to generate a 
vector of measured outputs from the array of photodetectors 
and to perform at least one of: 

analyzing properties of a target object; 

monitoring a change of the target object; and 

estimating a spectral profile of light from a target object 

through at least one of the following methods: 
linear, multi-linear, or nonlinear estimation, 
trained mapping; and 

pattern recognition. 

15. The method of claim 14, wherein the processing unit is 
configured to estimate a spectral profile of the light from the 
target object by sampling and digitizing a spectral response of 
the filters or combined spectral responses of each filter, and a 
spectral response of the photodetectors coupled to the filters, 
to thereby analyze properties of the target object. 


* * * * * 


